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Abstract 

o 

We propose a new type of gauge-invariant expansion of the ionization probability amplitudes of atoms by 
short pulses of electromagnetic radiation. Contrary to previous gauge-invariant approaches to this problem 
it does not require different partitions of the total Hamiltonian depending on the choice of gauge. In a 
natural way the atomic potential is treated as perturbation acting on an electron interacting with strong 
pulse. Whereas this is a standard assumption of strong field approximation (SFA), we show that grouping 
consequently together all terms of the same order in the atomic potential results in the expansion of the 
amplitude which is gauge invariant order by order, and not only in the limit of infinite series. In this approach, 
which is illustrated by numerical examples, the "direct ionization" and "rescattering" contributions are 
different from those commonly used in SFA - calculations. 

E 

1. Introduction 

Theoretical descriptions of the interaction of atomic and molecular systems with strong external fields 
is based mainly on the Keldysh-Faisal-Reiss (KFR) theory [TH3] where, in the first approximation, influ- 
ence of atomic forces on the electron dynamics is neglected, except of the exact initial state wave function. 
Assumption about dominance of the external field over atomic forces is justified in the case of very strong 
electric fields of the laser beam, basically of the order of the atomic electric field or higher. For a recent 
application of the Keldysh theory to the description of ionization by strong one-cycle pulses see, for in- 
stance, (3J. This type of approximation, known also as strong field approximation (SFA), has been applied 
to the description of interaction both with long optical wave-trains containing many cycles of the carrier 
wave [5J , and with relatively short few-cycle pulses [SJ [7] . Theoretical description in the framework of SFA 
including the rescattering contribution leads to a fair agreement with experimental data for above threshold 
ionization (ATI) spectra [9] . The KFR theory modified to account for asymptotic Coulomb effects leads 
also to a qualitative agreement with the measurements of momentum distributions in the direction of linear 
polarization of the pulse electric field [TOl [H] . 

One of main theoretical problems of SFA, despite of its well established status as a tool used for de- 
scription of interaction with strong fields, is connected with the proper choice of gauge of electromagnetic 
potential describing the field of the laser pulse. Two most popular gauges, the velocity gauge (VG) and 
the length gauge (LG), lead in general to different predictions [T2HT4] . In some cases the length-gauge 
SFA-calculation for detachment from negative ions matches better results obtained by numerical solution of 
time-dependent Schrodinger equation (TDSE) [T3], this gauge seems also to be "favored" by experimental 
results, giving predictions in agreement with measurements, at least quantitatively |15H17) . This conclusion 
has been questioned by Reiss [13] . who argued that velocity-gauge is more suitable for the description of 
experimental results for detachment fron negative fluorine ion, reported in |15j . Velocity gauge was earlier 
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advocated for by Cormier and Lambropoulos in |18j . The question of getting gauge invariant predictions in 
the case of interaction with few-cycle pulses in VG, LG and Henneberger frame [19] was discussed in the 
context of unitary transformations in [20j . 

Whereas gauge invariance of transition amplitudes (up to a phase factor) and probabilities does not raise 
questions as such, and can be easily proved formally, this fundamental property may be violated by various 
types of approximate calculations. In the theoretical description of the interaction with strong laser fields 
two approximations are used: dipole and SFA, and both can be considered as sources of gauge-dependence. 
Obviously, the question of "p-A" vs. "d-E" (VG vs. LG) description of interaction with electromagnetic field 
in dipole approximation is not new, especially in the context of bound-bound transitions in the perturbation 
theory (for reference, cf. eg. [2"TM27| ). In the case of bound- free transitions in strong laser fields it has been 
shown in |28j that VG - and LG - transition amplitudes are equivalent in all orders of the expansion provided 
appropriate initial and final state partitions of the Hamiltonian are chosen. In general, however, expansions 
of the amplitude in different gauges require various gauge dependent partitions of the Hamiltonian, and lead 
to the same result only in the limit of infinite series [29], and not "order by order" in the expansion. 

In this paper we propose an expansion of the amplitude which does not require any special, gauge 
dependent partitions of the Hamiltonian, and is manifestly gauge invariant in every order separately. The 
main idea is to collect together consequently all terms of a given order in atomic potential treated as 
perurbation in comparison with strong external field. The first gauge dependent term in standard SFA 
expansion is, in fact, of the first order in atomic potential, as can be seen for instance in |T3j [30] . We show 
in particular that next gauge dependent term in the expansion, corresponding to first order rescattering 
process [7], when combined with the first one, gives a complete gauge independent contribution of the first 
order in atomic potential, plus a term which is of next order in atomic potential. The iterative procedure 
leads then to an expansion which is manifestly gauge independent order by order. 



2. The manifestly gauge inavariant expansion of the amplitude 

The Hamiltonian describing interaction of an atomic system with external strong field in single active 
electron approximation has the following form (atomic units are used throughout) 

H(t) = H at + F(t), (1) 

where 

H a t = Y + V«ti ( 2 ) 

with V a t denoting time-independent atomic potential, F(t) describes the interaction with external field and 
p is the canonical momentum p = — iV. The Hamiltonian (JlJ can be equivalently written as 

H(t) = H F (t) + V aU (3) 

where Hp(t) is the Hamiltonian of an electron in the external field, ie. 

H F (t) = ^ + F(t). (4) 

Using partitions of the Hamiltonian defined by ([!]) and ^ one obtains two forms of the integral equation 
for the time evolution operator 

U(t,t') = U at (t,t')-i [ dtiUfatjF^UatihJ), (5a) 



U{t,t') = U F (t,t')-i / dhUfatJVatUp^t'), 



(5b) 



where U a t is generated by the atomic Hamiltonian ^ and Up - by the Hamiltonian Q. Transition amplitude 
from an initial bound electron state \<pi(t')) to a final continuum state \<fif(t)), orthogonal to the initial state 
(both are eigenstates of the atomic Hamiltonian), reads 

M={4 > f(t)\u{t,t l )\Ut')) 

dh{cj>f{t)\U{t,h)F{h)\Uh)), (6) 

where \<pi{t\)) = U a t(ti, t')\4>i(t')) . This exact expression for the amplitude is, up to a phase factor, inde- 
pendent on the choice of gauge. For short pulses the interaction Hamiltonian F(t) is practically zero for 
times earlier than some initial time U and later than a final time tf. Therefore, the time integration limits 
in ([6| can be replaced by U and t f . 

Expansion of the amplitude in powers of atomic potential can be achieved by consecutive iterations of 
(5b). First term corresponds to the replacement of the evolution operator by Up, so that 

M^M = -i / ' dti(tf(t)\U F (t,ti)F(h)\<lH(h)), (7) 



where also an approximation \(f>J(t)) for the final state has been used. Final state can be approximated by 
plane wave [31 [7], which is well justified for ionization of negative ions. In the case of ionization of a neutral 
atom one should use rather a plane wave distorted by asymptotic Coulomb phase [THl E] • The initial and 
approximated final state are not, in general, orthogonal to each other. Using equation fulfilled by Up, 



-id tl Up{t,t 1 ) = Up{t,t 1 ) 



r p 2 . 



(8) 



substituting the product UpF from Q into (J8|, performing integration by parts and using the equation 
fulfilled by \4>i{ti)), one finds 

M = (r f (t f )\U F (t f ,t t )\MU)) - (<t> a f (tf)\<t>i(tf)) 

dhi^itf^Upitf^VatlMh)}- (9) 



In the case of a periodic wave train or long pulses the boundary term gives a negligible contribution and the 
amplitude is then determined solely by the third term [31 H31 [3D] . 

Gauge transformation is implemented by a unitary operator of the form exp[i^ s (r, i)], where for trans- 
formations compatible with the dipole approximation Xg can be at most linear in the coordinate r. Vector 
potential, scalar potential and the evolution operator transform according to 

Ag = A-Vx g , (fig = ip + dtXg, 
U g (t,t') =e l ^^U(t,t')e- l ^ {t ' ) . (10) 

Transformed Hamiltonian has the form 

Hg{t) - e^^H^e-^^ - d tXg (r,t). (11) 



It follows from ( 11 ) that the interaction Hamiltonian in a new gauge reads 

F g (t) = F(t) - (l/2)(p • V X g + Vx s • P) 

+{l/2)(Vxg) 2 -d t Xg. (12) 

The Hamiltonian in a new gauge is again partitioned either according to ([!]) or according to ([3]), 



H g (t) = H at + F g (t) = H Fs (t) + V at . 
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(13) 



We shall consider a class of gauge transformations for which X g (t) — for t < ti and t > tf. This class of 
gauge transformations is characterized by 

X g (r,t)= 7 r- A(t), (14) 

where 7 is a real parameteiQ Whereas first two terms in (|9| are gauge-independent, the third one is not 
since U F (tf,t\) ^ U gF (tf,tij for transient times. Since the derivation leading from ([7| to ([9| can be done 
in any gauge, transition amplitude in the new gauge g can also be written as 



M g0 = (r f (t f )\u gF (t f ,u)\Mu)) - (r f (t f )Mtf)) 

dh^^UgFitf^VaMtl))- (15) 



u 



First two terms in (|9| and ( |15[ ) are equal but remaining terms differ among themselves, since for transient 
times U F {tf,t\) 7^ U gF (tt,ti). Next term, M\, in the expansion of the amplitude can be obtained from (5a| 



after substituting for U first iteration of ( 5b ) , which gives 

Ah = {-if f'dti (' dt 2 (r f {tf)\U F (t,h)V at 
Jti Jti 

xU F {t2,tx)F{ti)\4>ifa))' ( 16 ) 

Changing the order of integration, using again ^ to express the product UpF and performing integration 
by parts we obtain 



/ 



M x = i j dt 2 (r f (t f )\u F (t fl t 2 )v at \ut2)) 

dt 2 (^{t f )\U F {t f ,t 2 )V at U F {t 2 ,ti)\&{U)) + ---, (17) 

where omitted term is of the second order in atomic potential. Adding ^ and (17) cancels the gauge 
dependent contributions in both terms, leading to a manifestly gauge independent approximate expression 
for the amplitude of the form M w + where 

M<°> = (^(t f )\U F (t f ,ti)\MU)) - (<P}(tf)\Mtf)), (18a) 



M« = -i ^ dh^tfWFitfMWatUFit^U^U)). (18b) 
Jti 

3. Numerical calculations 

We have performed numerical calculations for a strong subfemtosecond laser pulse linearly polarized 
along the z axis. Vector potential in dipole approximation is 

TP 

A(t) = —ef(t) sm{ujt + ip), (19) 

UJ 

and E(t) = —d/dtA(t), Eq is the electric field amplitude, e is the polarization direction, w is the carrier 
wave frequency, ip is the carrier-envelope phase (CEP) and f(t) = exp [— i 2 /2r 2 ] is the envelope function. 
The value of r is 1.94 a.u., corresponding to FWHM of the intensity equal to 3.23 a.u. The carrier wave 



1 This is a subclass of a wider class of transformations considered in |29]. 
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length used in the calculations has been chosen as 72 nm (carrier wave period T = 2tt/oj = 240 as) and 
maximum electric field as Eg — 10 a.u., i.e. the maximum intensity I max ~ 7 x 10 18 W/cm . Possibilities 
of producing strong subfemtosecond pulses in a short wavelength region have been recently discussed eg. 
in [31H34]. 

We have calculated momentum distributions of final electrons for the detachment of H ion, modeled 
by a short range static potential of the form [7] 

V at {r) = -(d+ 9 -y-^ 1 (20) 

with g = d = 1 and /i = 1.56 a.u. This model potential supports a weakly bound state with the ionization 
energy 0.024 a.u., characteristic for the H ion. The initial bound state wave function has been found 
by numerical solution of the stationary Schrodinger equation and the amplitude was found by numerical 



calculation of the integrals ( 18 ) in momentum space. 

Results of calculation of the photoelectron spectra in the direction of polarization, together with the 
standard SFA results in the velocity gauge for comparision, are shown in Fig. [I] For a pulse with carrier- 
envelope phase (CEP) tp = (left panel) the spectra calculated with the use of M<°), and standard SFA - 
approach are symmetric around p z = 0, whereas the correction M' 1 ', to the amplitude, which includes part 
of a standard rescattering contribution, results in an asymmetric spectra. Spectra calculated with the use 
of the gauge-invariant approach (blue and red lines) differ substantially from the standard SFA results for 
smaller values of momentum. For tp = 7r/2 (right panel) the differences between those three approximations 
are much smaller. 



w(p z , 0) w(p z , 0) 




Figure 1: Normalized momentum distribution in the polarization direction as function of the longitudinal momentum for 
detachment from H - ion for ip = (left panel) and tp = tt/2 (right panel). Blue (dashed) line - distribution calculated using 
JW(°) (direct ionization), red (dotted) line - calculation with fJfW i.e with rescattering accounted for, green (continuous) 
line - standard SFA calculation in velocity gauge. 



Photoelectron spectra in the direction of polarization integrated over the transverse momentum, ie. 

1 f 00 

P(P*) = ^2 J Pt)Pt d Pt> ( 21 ) 

are shown in Fig. [2j Again, for <p = one can observe larger differences between gauge-invariant direct 
scattering, gauge-invariant rescattering and standard SFA approaches than in the case at tp = 7r/2, especially 
for smaller values of p z . For tp = 7r/2 the spectrum shows much richer structure than for tp = 0. As expected, 
for tp = 0, ie. even electric field, E(— t) — E(t), the "direct ionization" probablility cx |M'°)| 2 and the 
standard SFA probability are symmetric around p z = 0, as can be seen from left panels in Figs. [T] and [2j 
For tp = tt/2 and odd electric field this property is not true. 

Density plots of the probability distributions in the (p z , pt) plane are shown in Fig. [3] for the CEP 
tp = and in Fig. [4] for tp = n/2. In the first case one notes substantial differences between the distribution 
obtained from formla ( [18] ) and the standard SFA result, which could be also observed for this value of CEP 
in the p z distributions in the polarization direction. On the other hand, with CEP equal to 7r/2 results of 
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Figure 2: Normalized momentum distributions integrated over transverse momentum as function of longitudinal momentum for 
detachment from H - ion for ip = (left panel) and tp = 7r/2 (right panel). Blue (dashed) line corresponds to direct ionization, 
red (dotted) line takes rescattering contribution into account and green line represents conventional SFA in velocity gauge. 

the present calculations are very close to the standard SFA-approach, at least for the values of the pulse 
parameters used here. 








Figure 3: Normalized density plots of momentum distribution for detachment from H ion with tp = 0. Left panel - distribution 
calculated with the use of M(°) + MM, right panel - standard SFA approach 



G 




4. Conclusions 

We have obtained manifestly gauge invariant expressions for the amplitudes of ionization by short and 
intense laser pulses in the dipole approximation. The derivation is based on the SFA-approach but, contrary 
to its standard formulation, all terms of the same order in atomic potential were consequently grouped 
together. This procedure allowed to eliminate gauge-dependent terms resulting in a new type of expansion 
with redefined nontributions of the "direct ionization" and "rescattering" - type (M^ and M^ 1 ' respectively 



in (18 1). The approach proposed in this Letter is gauge-invariant order by order in the expansion for any 



gauge choice from the class of gauges defined by (14). This is in contrast with the approach used in 



for the standard SFA approach, where for a class of gauges compatible with dipole approximation gauge- 
invariance can be achieved in the limit of infinite series. The present approach is based on two "natural" 
partitions of the total Hamiltonian: atomic Hamiltonian plus interaction with the laser field (Eq. Q], or 
Hamiltonian of an electron in the laser field plus atomic potential [Eq. (|3j)]. These two partitions can be 



used for any gauge from the class (|l4|) with no necessity to use a notion of generalized field free Hamiltonian, 
as has been done in 



Numerical calculations of the detachment probabilities based on formulae ( 18 1 show for zero CEP sub- 
stantial differences in comparison with standard velocity gauge SFA, whereas for CEP equal to ir/2 the 
differences are much smaller, at leat for the pulse parameters used in the present calculations. The issue of 
a relation between the standard SFA/rescattering approach, as well as the problem of convergence, require 
further examination. Preliminary estimations suggest that present approach is more suitable for very short 
pulses, basically in the subfemtosecond region. 
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